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Abstract 

The scalar productions in heavy meson decays can provide a good platform to study not only heavy 
flavor physics but also their own physical properties in a dramatically different way. In this work, based on 
the assumption of two-quark structure of the scalars, the charmless hadronic B c — > SP, SV decays(here, 
S, P, and V denote the light scalar, pseudoscalar, and vector mesons, respectively) are investigated 
by employing the perturbative QCD(pQCD) factorization approach. In the standard model all these 
considered B c meson decays can only occur through the annihilation diagrams. From our numerical 
evaluations and phenomenological analysis, we find that (a) the pQCD predictions for the CP-averaged 
branching ratios(BRs) of the considered B c decays vary in the range of 10 -5 to 10~ 8 , which will be tested 
in the ongoing LHCb and forthcoming Super-B experiments, while the CP- violating asymmetries for these 
modes are absent naturally in the standard model because only one type tree operator is involved; (b) for 
B c — > SP, SV decays, the BRs of AS = processes are basically much larger than those of AS = 1 ones 
as generally expected because the different Cabibbo-Kobayashi-Maskawa(CKM) factors are involved; (c) 
analogous to B — > K*rf'> decays, Br{B c — > K + rf) ~ 5 x Br(B c — > n + n') in the pQCD approach, which 
can be understood by the constructive and destructive interference between the rj q and rj s contributions 
to the B c — > k + t] and B c — > k + t]' decays, however, Br(B c — » A"q(1430)^) is approximately equal to 
Br(B c —> Kq (1430)77') in both scenarios because the factorizable contributions from n s term play the 
dominant role in the considered two channels; (d) if ao(980) and k are the qq bound states, the pQCD 
predicted BRs for B c — > ao(980)(7r, p) and B c — > kK^*> decays will be in the range of 10~ 6 ~ 10~ 5 , which 
are within the reach of the LHCb experiments and could be measured in the near future; and (e) for the 
a (1450) and if,*(1430) channels, the BRs for B c a (1450)(7r, p) and B c if *(1430) modes in the 
pQCD approach are found to be (5 ~ 47) x 10~ 6 and (0.7 ~ 36) x 10~ 6 , respectively A measurement of 
them at the predicted level will favor the qq structure and help understand the physical properties of the 
scalars and the involved QCD dynamics in the modes, especially the reliability of the pQCD approach to 
these B c meson decays. 

PACS numbers: 13.25.Hw, 12.38. Bx, 14.40.Nd 
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I. INTRODUCTION 



The scalar mesons are especially important to understand because they have the same quan- 
tum numbers as the vacuum (J PC = 0++). Great efforts have been made by the physicists on 
both experimental and theoretical aspects to understand the inner structure of the scalars but 
it is well-known that the underlying structure of them is not yet well established (for a review, 
see e.g. [l|-0|)- Up to now, many different possible solutions to the scalars have been proposed 
such as qq, qqqq, meson-meson bound states or even supplemented with a scalar glueball. More 
likely, they are not made of one simple component but are the superpositions of these contents. 
The different scenarios tend to give very different predictions on the production and decay of 
the scalar mesons which are helpful to determine the dominant component. 

The first charmless B decay into a scalar meson, i.e., B — > fo(980)K, was measured by 
Belle jl] in 2002 (updated in [H]) and subsequently confirmed by BaBar 0] in 2004. After that 
these two B factories operated at KEK and SLAC respectively have found many decay channels 
with the scalars as one of the productions in B meson decays [H, 0] • These measurements should 
provide information on the nature of the scalar mesons. It is of enough reasons to believe that 
as a different unique insight to the internal structure of the scalars, the heavy B meson decaying 
into scalar mesons can provide a good place to explore their physical properties. 

Recently, the production of scalar mesons wit h q q structure in the two-body charmless B 
decays have been intensively studied in Refs. |8Hl2|] theoretically, in which many predictions 
are within the reach of the current B factory experiments and to be examined in the near 
future. It is hoped that through the study of B — )■ SP, SV(Here, S, P, and V are the light 
scalar, pseudoscalar, and vector mesons, respectively.) decays, old puzzles related to the internal 
structure and related parameters, e.g., the masses and widths, of light scalar mesons can receive 
new understanding. 

Experimentally, the Large Hadron Collider (LHC) experiment at CERN is running now, 
where the B c meson could be produced abundantly. Motivated by the forthcoming large number 
of B c production and decay events in the ongoing LHCb experiments, the scalar meson spectrum 
would become one of the most interesting topics for both experimental and theoretical studies 
in the near future. At that time, more and more channels with scalar mesons will be opened 
and got stringent tests from the experiments, which will help us to further explore the nature 
of the scalars. On the other hand, for B c meson, one can study the two heavy flavors b and c 
in a meson simultaneously. The B c meson decays may also provide windows for studying the 
perturbative and nonperturbative QCD, final state interactions, testing the predictions of the 
standard model(SM), and can shed light on new physics scenarios beyond the SM [131 ]. 

Inspired by the above observations, in this work, we therefore will focus on the two-body 
charmless hadronic decays B c — > SP, SV, which can only occur through the weak annihilation 
diagrams. The size of annihilation contributions is an important issue in B physics for many 
years. The importance of annihilation contributions has already been tested in the previous 



predictions of branching ratios of pure annihilation B — > D S K decays LJ], direct CP asymme- 



tries of B° — y 7r + 7r , K + tt decays 15l4l7l] and in the explanation of B — > <pK* polarization 



problem [la . Il9 ] though there still exist much different viewpoints 1 . The two-body B decays 



1 Recently, the authors announced in Ref. [2p| ] that the annihilation contributions in charmless hadronic B 
decays are real and small in the soft-collinear effective theory [2l[ at leading power. While the authors in 
another work [22l | discussed that they may be the almost imaginary contributions, which can generate a sizable 
strong phase. This discrepancy between these two approaches/methods needs to be clarified definitely by the 
experiments in the future. 
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into the final states with one scalar meson may suffer from large weak annihilation contributions, 
which have been analyzed in Refs. 10h12|] preliminary. Thus it is very interesting to explore 
the size of annihilation contributions in these considered B c — > SP, SV channels, which will also 
be helpful to investigate the annihilated decay mechanism and the physical properties of the 
scalars. 

In this paper, we will study the CP-averaged branching ratios(BRs) of charmless hadronic 
B c — > SP, SV decays by employing the low e nergy effective Hamiltonian [23| and the perturba- 
tive QCD(pQCD) factorization approach 15], I la. |24|]. By keeping the transverse momentum kx 
of the quarks, the pQCD approach is free of endpoint singularity and the Sudakov formalism 
makes it more self-consistent. Rather different from the QCD factorization approach [25] and 
soft-collinear effective theory, the pQCD approach can be used to calculate the annihilation 
diagrams straightforwardly 26)], as have been done for example in Refs. IHHElHIzJSl. 

The paper is organized as follows. In Sec. UU we present a brief review of light scalar mesons 
and the formalism of pQCD approach. The wave functions and distribution amplitudes for heavy 
B c and light scalar, pseudoscalar, and vector mesons are also given here. Then we perform the 
perturbative calculations for the considered B c — > SP, SV decay channels with pQCD approach 
in Sec. IIIII The analytic formulas of the decay amplitudes for all the considered modes are 
also collected in this section. The numerical results and phenomenological analysis are given in 
Sec. IIVI Finally, Sec. [V] contains the main conclusions and a short summary. 



II. LIGHT SCALAR MESONS, FORMALISM, AND WAVE FUNCTIONS 



A. Light scalar mesons 



Up to now, the people have discovered many scalar states experimentally but know little 
about their underlying structures, which are not well established theoretically yet (for a review, 
sec Refs. UJH). 

According to the meson particle collected by the Particle Data Group [lj], the 
light scalar mesons below or near 1 GeV, including ao(980), l£g(800)(or k), /n (600) (or a), and 
/o(980), are usually viewed to form an SU(3) flavor nonet; while scalar mesons around 1.5 GeV, 
including o (1450), #5(1430), / (1370), and / (1500)// (1710), form another nonet 2 . 

Recently, Cheng, Chua, and Yang [10j proposed two possible scenarios to describe these light 
scalar mesons in the QCD sum rule method: 

1. In scenario 1(S1), the scalar mesons in the former nonet are treated as the lowest lying 
states, and in the latter one as the corresponding first excited states, respectively. Based 
on the naive two-quark model, the flavor structure of the light scalar mesons in SI read 

a = —= (uu + dd) , fo = ss, 
V2 

clq = ud , ag = — t= (uu — dd) , clq = du , 
v2 

k + = us , k° = ds , R° = sd , K~ = su . (1) 



2 For the sake of simplicity, we will use ao and fo to denote ao(980) and /o(980), respectively, unless otherwise 
stated. We will also adopt the forms a, Kq, f, and /' to denote the scalar mesons ao(980) and ao(1450), Kg (800) 
and 7^(1430), / (600) and / (1370), and / (980) and / O (1500)// O (1710) correspondingly in the following 
sections, unless otherwise stated. 
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Here, it is assumed that the lightest a and heaviest fo in the lighter scalar nonet has the 
ideal mixing. But various experimental data indicate that fo should not have the pure 



ss component and the isoscalars a and fo must have a mixing of /q and /q 10J, which is 
analogous to ij — 7/ mixing system, 

fo ) ~ \ sin 6» cos O J \ / S 
with /q = (uu + dd)/V2 and /q = ss, where 0$ is the mixing angle between a and fo. 



Many works have been made to explore the mixing angle 6q [31[: #0 hes in the ranges of 
25° < 6»o < 40° and 140° < 6> < 165°. But the fact that <9 tends to be not a unique 
value, which indicates that a and fo may not be purely qq states. 

While for the mixing of the isosinglet scalar mesons .fn(1370), fo(1500), and /o(1710), 
which have been discussed in detail in the literatures (See Ref. [32| and references therein). 



In this work, we will adopt the mixing mechanism as given in Ref. 32] 



/o(1370) \ 

/o(1500) = ( -0.54 0.84 0.03 | ( / S | . (3) 
/o(1710) / 




As discussed in [32|, it is evident that /o(1370) and /o(1500) mainly consists of /q and /5 s , 
just with small or tiny glueball components, however, /o(1710) is composed primarily of 
the scalar glueball, i.e., fff. We will therefore only take the scalar mesons /o(1370) and 
/o(1500) into account in the present work, and leave the contribution from scalar glueball 
content for future study. 

2. In scenario 2(S2) that the scalar mesons in the latter nonet are the lowest lying resonances 
and the corresponding first excited states lie between (2.0 ~ 2.3) GeV. S2 corresponds to 
the case that light scalar mesons below or near 1 GeV are four-quark bound states, while 
all scalar mesons are made of two quarks in SI. In order to give quantitative predictions, 
since we do not know how to deal with the four-quark states in the factorization approach 
presently, we here just consider the evaluations on the scalar mesons with qq structure in 
S2. 

In short, we will investigate these light scalar mesons in the pure annihilation B c — > SP, SV 
decays with the assumption of two-quark structure proposed in the above two possible scenarios. 



B. Formalism of pQCD approach 

Since the b quark is rather heavy, we work in the frame with the B c meson at rest, i.e., with 
the B c meson momentum P\ = ^=-(1,1, Ot) in the light-cone coordinates. For the charmless 

hadronic B c — > M2M3 3 decays, assuming that the M2 (M3) meson moves in the plus (minus) 
z direction carrying the momentum P2 (P3) and the longitudinal polarization vector (eg)(if 
M 2 (3) is the vector meson). Then the two final state meson momenta can be written as 

"•(1 /i/iO/). P 3 = ^§( r |,l- r 2,0T), (4) 



3 For the sake of simplicity, we will use M2 and M3 to denote the two final state light mesons respectively, unless 
otherwise stated. 
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respectively, where r 2 = mM 2 / rn B c an d r 3 = mM 3 /mB c - When M2 or M 3 is a vector meson, the 
longitudinal polarization vector, e^, or e 3 j can De given by 

= ^Ml-ri,-r 2 2 ,0 T ), or 4 = -^(-r 3 2 , 1 - r 2 2 , T ). (5) 
V 2mM 2 V2mM 3 

Putting the (light-) quark momenta in B c , M2 and M3 mesons as k\, fc 2 , and ^3, respectively, 
we can choose 

fci = (a; 1 i^-,0,kir) J k 2 = (x 2 P+, 0, k 2T ), k 3 = (0, x 3 P~, k 3T ). (6) 

Then, for B c — > M2M3 decays, the integration over fcf, k% , and A;^ will conceptually lead to 
the decay amplitudes in the pQCD approach, 



A(B C — > M 2 M 3 ) ~ Jdxidx2dx 3 bidbib2db2b 3 db 3 



•Tr 



C{t)<S> Bc {x 1 ,b 1 )<S>M 2 (x2,b2)<S>M 3 {x3,b 3 )H(x l ,b l ,t)S t (x i )e- s{t) ,(7) 



where bi is the conjugate space coordinate of kiT, and t is the largest energy scale in function 
H(xi,bi,t). The large logarithms ln(mw/t) are included in the Wilson coefficients C{t). The 
large double logarithms (ln 2 Xj) are summed by the threshold resummation [3^], and they lead 
to St(xi) which smears the end-point singularities on X{. The last term, e~ s ^\ is the Sudakov 



form factor which suppresses the soft dynamics effectively 34]. Thus it makes the perturbative 
calculation of the hard part H applicable at intermediate scale, i.e., niB c scale. We will calcu- 
late analytically the function H(xi,bi,t) for the considered decays at leading order(LO) in a s 
expansion and give the convoluted amplitudes in next section. 

For these considered decays, the related weak effective Hamiltonian H e g [23| can be written 

as 

H cS = ^ [V* b V uD (C 1 ( M )0 1 ( M ) + C 2 (M)0 2 (/i))] , (8) 

with the local four-quark tree operators Oi ; 2 

01 = %7 M (1 - 75) Dcfip^il - j 5 )b a , 

02 = %7 M (1 - J 5 )Dpc a ^(l - J5 )b a , (9) 

where V c b,V u o are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, "I?" denotes the 
light down quark d or s and Cj(//) are Wilson coefficients at the renormalization scale /x. For 
the Wilson coefficients Ci^il^), we will also use the leading order expressions, although the 



next-to-leading order calculations already exist in the literature 23J]. This is the consistent way 
to cancel the explicit /i dependence in the theoretical formulae. For the renormalization group 
evolution of the Wilson coefficients from higher scale to lower scale, we use the formulas as given 
in Ref. [3] directly. 



C. Wave functions and distribution amplitudes 

In order to calculate the decay amplitude, we should choose the proper wave function of 
the heavy B c meson. In principle there are two Lorentz structures in the B q (q = u, d, s) or B c 
meson wave function. One should consider both of them in calculations. However, since the 
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contribution induced by one Lorentz structure is numerically small [28|, [35] and can be neglected 



approximately, we only consider the contribution from the first Lorentz structure. 

<D Bc (z) = -^=[(r+ MbJ^b^x)}^ . (10) 

Since B c meson consists of two heavy quarks and niB c — nij, + m c , the distribution amplitude 
§B C would be close to 5{x — m c /rriB c ) in the non-relativistic limit. We therefore adopt the 
non-relativistic approximation form of 4>b c as [361 ] . 

<t>B c (x) = ^==S(x - m c /m B J , (11) 

where fs c and N c are the decay constant of B c meson and the color number, respectively. 
The wave function for the scalar meson(S') can generally be defined as, 

d>s(*) = -^=Jr<Ps(x)+rns4(x)+m s (^-l)4(x)} ap , (12) 

where (ps and 4> s ' , and mj are the leading twist and twist-3 distribution amplitudes, and mass 
of the scalar meson, respectively, while x denotes the momentum fraction carried by quark in 
the meson, and n = (1, 0, Ot) and v = (0, 1, Ot) are dimensionless light-like unit vectors. 

In general, the leading twist light-cone distribution amplitude </>s(x,/i) can be expanded as 
the Gegenbauer polynomials 1(], 37] : 



0s(x, M ) = -^x{l-x){f s ^) + h^)^BMCU 2 {2x-l)Y (13) 

where fsW and fs(l^), B m (fi), and C m (t) are the vector and scalar decay constants, Gegen- 
bauer moments, and Gegenbauer polynomials for the scalars, respectively. 

Because of the charge conjugation invariance, neutral scalar mesons cannot be produced by 
the vector current and thus 

U = ffo = f a = 0. (14) 

For other scalar mesons, there exists a relation between the vector and scalar decay constants, 

fs = Vsfs and fi s = -. -r-- , (15) 

where mi and rri2 are the running current quark masses in the scalars. For the neutral scalar 
mesons /o, Q-o and a, fs vanishes, but the quantity fs = fsHs remains finite. 

The values for scalar decay constants and Gegenbauer moments in the scalar meson distri- 
bution amplitudes have been investigated at scale fi = 1 GeV in Ref. [l( 



f ao = 0.365 ± 0.020 GeV, B x = -0.93 ± 0.10, B 3 = 0.14 ± 0.08 (SI), 
f R = 0.340 ± 0.020 GeV, B x = -0.92 ± 0.11, B 3 = 0.15 ± 0.09 (SI), (16) 
//„ = 0.370 ± 0.020 GeV, B\ = -0.92 ± 0.11, B 3 = 0.15 ±0.09 (SI); 

/ O0 (l450) = -0.280 ± 0.030 GeV, B x = 0.89 ± 0.20, B 3 = -1.38 ± 0.18 (SI) , 
/a (1450) = °- 460 ± °- 050 GeV > B i = -°- 58 ± °- 12 > ^3 = -°- 49 ± °- 15 (S2) ; (17) 
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FIG. 1: Typical Feynman diagrams for two-body charmless hadronic B c — > SP(PS) decays at leading 
order. By replacing the pseudoscalar meson P with the vector meson V, which will lead to the diagrams 
for B c -> SV(VS) modes. 

7^(1430) = -0.300 ± 0.030 GeV, B\ = 0.58 ± 0.07, B 3 = -1.20 ± 0.08 (SI) , 
/k*(1430) = 0.445 ± 0.050 GeV, B x = -0.57 ± 0.13, B 3 = -0.42 ± 0.22 (S2) ; (18) 



// (1500) = -0.255 ± 0.030 GeV, B x = 0.80 ± 0.40, B 3 = -1.32 ± 0.14 (SI) , 
// (i500) = 0.490 ± 0.050 GeV, B x = -0.48 ± 0.11, B 3 = -0.37 ± 0.20 (S2) . (19) 

As for the twist-3 distribution amplitudes 4>g and (fig, we adopt the asymptotic forms: 

* i = vk /s ' * l= ijm h[l - 2x) - (20) 

Note that for the distribution amplitudes of strange scalar meson, x stands for the momentum 
fraction carrying by s quark. 

For pseudoscalar meson(P), the wave function can be generally defined as, 

$ P (x) = -^=75 {ft>p(x) + m^P(x) + m^(# - l)^(x)}^ , (21) 

APT p 

where (f> p ' ' and tuq are the distribution amplitudes and chiral scale parameter of the pseu- 
doscalar meson, respectively. 

For the wave functions of vector meson(F), one longitudinal(L) polarization is involved, and 
can be written as, 

where eh denotes the longitudinal polarization vector of vector mesons, satisfying P ■ e = 0, 

t s 

4>v and <j)y , and my are the leading twist and twist-3 distribution amplitudes, and mass of 

APT 

the vector meson, respectively. For the distribution amplitudes of pseudoscalar <j) p ' ' , and 
longitudinal polarization, (f>y and (py S to be used in this work, we will adopt the same forms as 
that in the literatures (See Ref. |30(j and references therein). 



III. PERTURBATIVE CALCULATIONS IN THE PQCD APPROACH 

From the effective Hamiltonian ([8]), there are 4 types of diagrams contributing to the B c — > 
M2M3 decays as illustrated in Fig. [TJ which result in the Feynman decay amplitudes j^ 2Ma 
and A4na M:i i where the subscripts fa and na are the abbreviations of factorizable and non- 
factorizable annihilation contributions, respectively. Operators 0\^. are (V — A)(V — A) currents, 
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we therefore can combine all contributions from these diagrams and obtain the total decay 
amplitude as, 

A(B C -> M 2 M 3 ) = V: b V uD {fB c Ff a 2Kh ai + M% Ma a} , (23) 

where a\ = Ci/3 + C 2 . In the next two subsections we will give the explicit expressions of 
jtM 2 M 3 ^ M M2M Z anc j t h e decay amplitude A(B C -> M 2 M 3 ) for 5 C M 2 M 3 decays: including 
32 £ c -> SP(PS) and 30 £ c -> Sy(FS) decay modes. 



A. B c -> SP(PS) decays 

In this subsection, we will present the factorization formulas for charmless hadronic B c — > 
SP(PS) decays. From the first two diagrams of Fig. [TJ i.e., (a) and (b), by perturbative QCD 
calculations, we obtain the decay amplitude for factorizable annihilation contributions as follows, 

rl poo 

Tfa = 8nCFm Bc / dx 2 <ix 3 I b2db2b 3 db 3 
c Jo Jo 

X {hf a (l ~ x 3, %2, h, b 2 )Ef a (t a ) [x 2 (f>s(x 2 )(l)p(x 3 ) + 2r SVq '<t>p(x 3 ) 

x ((x 2 + l)0|(x 2 ) + (a?a " 1)05(^2))] + h fa (x 2 , 1 - x 3 , b 2 , h)E fa (t b ) 

x [(x 3 - l)^ s {x 2 )^{x 3 ) + 2r s r£ 4{x 2 ) {(x 3 - 2)4> p P {x 3 ) - x 3 cf>l(x 3 ))] } , (24) 

where 4>s(p) corresponds to the distribution amplitudes of mesons S(P), rs = ms/rriB c , Tq = 
m[/mg c , and Cp = 4/3 is a color factor. The function hf a , the scales ti and Et a (t) can be 



found in Appendix B of Ref. 30(] . 



For the nonfactorizable diagrams (c) and (d) in Fig. [H all three meson wave functions are 
involved. The integration of b 3 can be performed using 5 function 5(b 3 — 6 2 ), leaving only 
integration of b\ and 6 2 . The corresponding decay amplitude is 



■Mna = —^—^Cprn Bc I dx 2 dx 3 f bidbib 2 db2 

Jo Jo 

x {Ka( x 2, x 3 ,b 1 ,b 2 )E na (t c ) [(r c - x 3 + l)(f)s(x 2 )(f>p(x 3 ) + r s r p (0f(x 2 ) 
x ((3r c + x 2 - x 3 + l)4>p(x 3 ) - (r c - x 2 - x 3 + l)<fp{x 3 )) + <^s(x 2 ) 
x ((r c - x 2 - x 3 + l)0£(x 3 ) + (r c - x 2 + x 3 - l)0^(x 3 )))] - E na (t d ) 
x [(r 6 + r c + x 2 - l)<fe(x 2 )0p(x 3 ) + r s r^ (0f (x 2 )((4r b + r c + x 2 - x 3 
-1)</>p(x 3 ) - (r c + x 2 + x 3 - l)4> T P {x 3 )) + 4> T s {x2){{r c + x 2 + x 3 - 1) 
x^(x 3 ) - (r c + x 2 -x 3 - l)4>p(x 3 )))} hi a (x 2 ,x 3 ,biM)) , (25) 

where r b = mb/mB c , r c = m c /mp, c , and r b + r c ~ 1 in l? c meson. 

Likewise, we can get the analytic factorization formulas of the contributions from l? c — > PS 
decays easily, 

pi poo 

Ffa = 87rCpm| c J dx 2 dx 3 J b 2 db 2 b 3 db 3 

x {hf a (l - x 3 ,X2,b 3 ,b 2 )Ef a (t a ) [x2<pp{x2)(t)s{x 3 ) - 2r^rs^f(x 3 ) 

x ((x 2 + i)<ffp(x 2 ) + (x 2 - l)0p(x 2 ))] + h fa (x 2 , 1 - x 3 , b 2 , b 3 )E fa (t b ) 

x [(x 3 - l)<f4(x 2 )M^) ~ 2r p r 5 ^(x 2 ) ((x 3 - 2)0f (x 3 ) - x 3 <j%(x 3 ))] } , (26) 
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■Mna = —ir^ 7r CFrn 2 Bc f dx 2 dx 3 I b\db\b 2 db 2 

Jo Jo 

x {h c na (x 2 ,x 3 ,b 1 ,b 2 )E na (t c ) [(r c - x 3 + 1)<Pp(x 2 )4>s{x 3 ) - r^r s {4>p(x 2 ) 

x ((3r c + x 2 - x 3 + l)#f (x 3 ) - (r c - x 2 - x 3 + 1)</>s(x 3 )) + 4>l(x 2 ) 
x ((r c - x 2 - x 3 + (x 3 ) + (r c - x 2 + x 3 - l)0p(x 3 )))] - E na (t d ) 
x [(r 6 + r c + x 2 - l)(f>p(x 2 )<f>s{x3) ~ r^rs ((j)p(x 2 )((Ar b + r c + x 2 - x 3 
-l)(f>S (x 3 ) - (r c + x 2 + x 3 - \)4> T s {x 3 )) + (j)p(x 2 )((r c + x 2 + x 3 - 1) 
x0|(x 3 )-(r c + X2-x 3 -l)0s(a;3)))]^ a (x 2 ,X3,6i,6 2 )} . (27) 

Based on Eqs. (|23M27]h we can write down the total decay amplitudes for 32 B c — > SP(PS) 
decays straightforwardly, 

A(B c ^a+7r°) = V: b V ud {[f Bc F]la 1+ MZCi] 

-lfB c ^ ai + MZCi]}/V2, (28) 
A(B c ^a\+) = ^{[/b/J>i+^K] 

-[fB c Tf>x+MZCi]}/V2; (29) 

A(B C -> a+r?) = V^Kd {[/b c ^ 9 «i + AOd] 

+[/B c ^7r«i + Ml q :d] } cos , (30) 

+ [fB c j rV fa a ai + M^aCi] \ sin ; (31) 



„4(5 C -> /tt+) = V* b V u d | [fBcFj^ai + Ml f Jd] 

+ [fs^fa^ + } cos O , (32) 

+ [fs^ax + MStc x ] } sin 6 ; (33) 

„4(5 C -»• K* + if°) = V^Kd {/b^^i + M^FCi} , (34) 

^(B c ^V+) = VSKd {/b c -F5' A ' + «i + Mfs K+ C x ) ; (35) 

,4(B C -> K*°7r+) = KfcK, {/b^V + A^Ci} , (36) 

= V2.4(5 C -> ^ * + 7T ) ; (37) 

A(B C -> a + K°) = V: b V us {f Bc Ff a a ai + M%&} , (38) 

= V2A(B C -»• i-T+a ) ; (39) 
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A{B C K*+rj) 



{/b c 
V4K S {/b c 

+ 



T' a "sin, 



■M^°sm< 



C'i 



}■ 



Ff° Vq sin ^ + J r J a K *° cos. 



Oi 



A^na 077 * Sin + cos < 



Ci} ; 



(40) 
(41) 



,4(B C -> /if H 



A(B C -> f'K + ) 



cos O 



A^n/° cos 6» - X{° a ^ sin O 



-^/f sin O 



+ 



■}• 



"^ /o %m0o + ^cos0o 



+ 



Mna^ Sin O + Mn°7 COS O ] d } 



ai 



(42) 



(43) 



B. £ c -> 51/(^5) decays 

After the replacement of the pseudoscalar meson P with the vector meson V in Figure [H 
we will get the Feynman diagrams for pure annihilation B c — > SV(VS) modes at leading order. 
By following the same procedure as stated in the above subsection, we can obtain the analytic 
decay amplitudes for B c — > SV decays, 



-pSV 



M 



SV 



pi poo 

—8irCF'm'^ c / dx 2 dx 3 j b 2 db 2 b 3 db 3 
c Jo Jo 

x {hfai} ~ x 3 ,x 2 ,b 3 ,b 2 )E fa (t a ) [x 2 <ps(x2)4>v{x3) - 2r s r V (/) s v (x 3 ) 

x ((x 2 + l)0|(x 2 ) + (X2 - l)(f) T s{x 2 ))} + h fa (x 2 , 1 - x 3 , b 2 , b 3 )E fa (t b ) 

x [(x 3 - l)(f> s (x 2 )<f> v (x 3 ) - 2r s r V (f> s s (x 2 ) ((x 3 - 2)<j> s v (x 3 ) - x 3 <p v {x 3 ))] } , (44) 

16\/6 f 1 f°° 

— irCFm 2 Bc I dx 2 dx 3 I b\db\b 2 db 2 

3 c Jo Jo 

x {Ka(x2,x 3 ,b 1 ,b 2 )E na (t c ) [(r c - x 3 + l)4>s(x 2 )4>v(x 3 ) - r s r v (</>f (x 2 ) 
x((3r c + x 2 - x 3 + l)0y(x 3 ) - (r c - x 2 - x 3 + l)(j>y{x 3 )) + (f>s(x 2 ) 
x((r c - x 2 - x 3 + 1)^(^3) + (r- c - x 2 + x 3 - l)^(x 3 )))] - £ n a(*d) 
x [(r b + r c + x 2 - 1)05(^2)^(^3) - rsry (</>f (x 2 )((4r 6 + r c + x 2 - x 3 
-l)<p s v (x 3 ) - (r c + x 2 + x 3 - l)<^(x 3 )) + 4>s(x 2 )((r c + x 2 + x 3 - 1) 
^^v(x3)-(rc + x 2 -x 3 -l)^ t v (x 3 )))]h^ a (x 2 ,x 3 ,b 1 ,b 2 )^ , (45) 

with ry = my /wib c - 

Similarly, the factorization formulas for B c — > VS decays can be easily obtained but with the 
simple replacements in Eqs. (|44|45p as follows, 



4>s < — > <Pv, 



< — > 



< — ► 4> v , r s i — >■ r v 



(46) 
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The total decay amplitudes of the 30 B c — > SV(VS) decays can therefore be written 

A(B c ^a + p°) = v: b V ud {[f Bc T a / a a 1+ MZCi} 
-[fB^ar + MZCi]}/^, 

A(B c ^a°p + ) = V: b V ud {[f Bc T p f a a a 1+ MZCi] 
-[fB^f a ai+MZCi]}/V2; 



A{B c ^a+u) = V: b V ud {[f Bc T^a 1+ MZCi] 
+[f Bc jySa 1+ MZCi]}/V2; 



A(B C -> 



A(B C -> f' P +) 



v: b V ud { [f Bc T p /J ai + MtSd] 
+ [/B c ^ P ai + Mttd}} cos0 o , 
v: b Vud { [fB c T p /J ai + M^d] 
+ [f Bc T f f i P ai + M^fCi]} sin d ; 



A(B C ^K*+K*°) = V; b V ud U Bc F^a 1 + Mna^CA , 
A(B c ^Kl°K* + ) = V: b V ud {f Ba jf} K *\ 1 +MB K * + C 1 } ; 



K K * 



A(B c ^K*°p+) = V; b V us {f Bc Tf a 5p a 1 + Mni P C 1 ] , 
= V2A(B C -> X* V) ; 

«4(-Bc — > a + K*°) = V: b V us {f Bc Tf; a a 1 +M^ a C 1 } , 
= V2A{B C ^ K* + a°) ; 



-K*L 



A(B C -> = F C * 6 F US (/s c ^; a ai + MnaCij /V2 , 



^(B c -> f'K* 



FfJ" cosOo-Fjf sin O 



+ [Mnl fS cos O - sin O 



C 



sm9 + F f ° a cos O 



+ [mTJ" sin O + 7W£ X * cos ( 



ai 



ai 



C; 



}■ 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



In this section, we will make the theoretical predictions on the CP-averaged BRs for those 
considered B c — > SP, SV decay modes. First of all, the central values of the input parameters 
to be used are given in the following, 

• Masses (GeV): 

niB c = 6.286 , mj, = 4.8 , m c = 1.5 ; 

m K * = 0.892 , m p = 0.770 , m w = 0.782 ; 

m K = 0.800, TTi/(j —— 0.600 , m /o = 0.980; 
m a (i450) = 1-474 , m K .( 1430 ) = 1.425 , m /o{1370) = 1.350 , m /o(1500) = 1.505 ; 

m$ = 1.6 , m v Q q = 1.08 , m v s = 1.92 . (62) 

Decay constants (GeV): 



mw 


= 80.41 


771,0 


= 1.02 , 


m ao 


= 0.985 


a (1450) 


= 1.474 


< 


= 1-4, 



U = 0.231 , fj = 0.200 , f K * = 0.217 , f K , = 0.185 ; 
f p = 0.209 , /J = 0.165 , f u = 0.195 , /J = 0.145 ; 
U = 0-131 , f K = 0.16 , f Bc = 0.489 . (63) 

• QCD scale and B c meson lifetime: 

A^T 4) = 0.250 GeV, t Bc = 0.46 ps. (64) 

Here, we adopt the Wolfenstein parametrization, and the updated parameters A = 0.814, 
A = 0.2257, p = 0.135, and f) = 0.349 \1\ for the CKM matrix. In numerical calculations, central 
values of the input parameters will be used implicitly unless otherwise stated. 

For B c — > SP, SV decays, the decay rate can be written as 

T = " r 2 s )\A(B c -> M 2 M 3 )| 2 , (65) 

where the corresponding decay amplitudes A have been given explicitly in Eqs. (|28H43p and 
Eqs. (|47H6ip . Using the decay amplitudes obtained in last section, it is straightforward to cal- 
culate the CP-averaged BRs with uncertainties as presented in Tables Ul lVlIII The dominant 
errors come from the uncertainties of charm quark mass m c = 1.5 =t 0.15 GeV, the scalar decay 
constants fs, the Gegenbauer moments en of the relevant pseudoscalar or vector meson distribu- 
tion amplitudes, the Gegenbauer moments B{ of the scalar meson distribution amplitudes, and 
the chiral enhancement factors TOq = 1.4 ± 0.3 GeV and uiq = 1.6 ± 0.1 GeV, respectively. 

Among the considered B c — > SP, SV decays, the pQCD predictions for the CP-averaged BRs 
of those AS = processes are basically much larger than those of AS = 1 channels (one of the 
two final state mesons is a strange one) , the main reason is the enhancement of the large CKM 
factor iKtd/Kjsl 2 ~ 19 for those A5 = decays as generally expected. Maybe there exist no 
such large differences for certain decays, which is just because the enhancement arising from the 
CKM factor is partially cancelled by the difference between the magnitude of individual decay 
amplitude. The pQCD predictions for the CP-averaged BRs of considered B c decays vary in 
the range of 10 -5 to 10 -8 . For B c — > ao(1450) + 7r° decay with rate of 10~ 5 ~ 10~ 6 for example, 
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TABLE I: The pQCD predictions of branching ratios(BRs) for the AS = processes of charmless hadronic 
B c — >• (clq, k, a, /o)(tt, K, rj, rj') decays in SI. The source of the dominant errors is explained in the text. 



AS = 
Decay modes 


BRs (10~ 6 ) 


B c -> a+7T° 
B c -> ag7r+ 


c c+2.3/„ \+0.9/? \+2.1/_7r\+1.4/DS , \+0.4/ \ 
b.5_! . 5 (mcJ_o.6l7Sj-1.4(«2j-l.l(- t? l,3j-0.7l m Oj 
o c+1.6/_ \+0.4/? \ + 1.0/ ttN+I.I/rS \+0-7(V„ \ 


B c -> 77 x 10 
B c -> a^' x 10 


3.6l^(m c )l°;|(/5)lS;I(^)lJ; 6 6 (^ 3 )±H(^o) 
2.4^(m c )t°; 3 2 (/ s )iS; 3 (ar)±Si(Sf.3)^(mo) 


B c -> k°K+ 
B c ->■ K°k+ 


q /1+2-V™ A+°- 5 ('f \+l-8/'„K \ + 14/nS \+0.2/_, \ 

2.1^(m c )t^(/5)iS; 5 1 (af 2 )ig;I(Sf. 3 )i^(mo) 


B c -> 7T+CT x 10 
B c -> 7T+fo x 10 


3.2^K)t°i(/ s )il;3(aI)l°; 9 7(Sf i 3)^: 3 (rr l o)(fo q ) 

i.8tUK)^:2(/s)i^K)^i(sf, 3 )^(rno)(fo q ) 



we show the decay amplitudes arising from both factorization and nonfactorization annihilation 
contributions explicitly(in unit of 10~ 3 GeV 3 ), 

A fa (B c -»■ a (1450) + 7r°) = 0.292 + i2.489; A na (B c -> a (1450) + 7r°) = 6.717 + i7.508(66) 

in SI, while 

A fa (B c -> a (1450) + 7T ) = 0.553 - i0.356; A na {B c -»■ a (1450) + 7r°) = 3.161 - »5.137(67) 

in S2, where the central values are quoted for clarification. One can find that the dominant non- 
factorizable decay amplitude governs this channel and subsequently results in the large branching 
ratio in both scenarios, which can be seen in Table HTT1 The other modes with large decay rates 
can be analyzed similarly. 

As discussed in Ref. [38J, the B c decays with the branching ratio of 10~ 6 can be measured at 
the LHC. Hence our pQCD predicted BRs with 10~ 6 or larger for these B c — > SP, SV decays are 
expected to be measured in the ongoing LHCb experiments, which will be very helpful to study 
the physical contents of the scalars and the involved QCD dynamics and annihilation mechanism 
in the these considered channels. Moreover, there is no CP violation for all these decays within 
the SM, since there is only one kind of tree operator involved in the decay amplitude of all 
considered B c decays, which can be seen from Eq. (|23p . 



A. B c -> a (P, V) and B c -> a (1450)(F, V) decays 

In this subsection, we will make some discussions on the B c — > a(P, V) decays involving 14 
AS = and 8 AS = 1 processes, respectively. 

From the numerical results for considered modes as given in the Tables HI IIII|, IV l and I VIII one 
can find that the CP-averaged BRs for all the AS = B c — > a(P, V) processes are in the range 
of 10 -6 ~ 10 -5 within the theoretical errors except for B c — V a^rf^ decays, which are expected 
to be tested by the ongoing LHCb measurements and the forthcoming Super-B experiments. 
Since we make the perturbative calculations based on the assumption of two-quark structure for 
the scalars, once these theoretical predictions could be verified by the related experiments, then 
these results will help us to explore the underlying structure of the scalar a meson. 
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TABLE II: Same as Table U but for the AS = 1 processes of charmless hadronic B c — > 
(a , k, a, f )(ir, K,rj,r}') decays in SI. 



A5 = 1 
Decay modes 
B c -+ a+K° 



a°IC 



BRs (10~ 7 ) 

(m c )tg: 



4.0 
2.0 



+0.4 
0.9 

+0. 




2(me)±S: 



0.7/nS ^1+0-0/ \ 
1.0l- tf l,3j-0.2l m 0j 

2/f \+0.9/„K ^+0.4/ R S \+0.0/„, \ 
3USj_0.9l a l,2j-0.5l ±f l,3J-0.ll TO 0j 



B c 



B c -> k+t?' x 



10 



4.5 
8.8 



3/„. \+0.l 
9 (m c )7- 



+2.4, 



-0 

4/ \+1.3/ J \+0.6/ n' 



(/s)tS 



Ml 



0.9 



-0.0/ \ 

-o.o(™o) 



+3.7/ R S \+0.0/ \ 



B c ->■ K°7T+ 
B r -> K+7T 



z - 1 -o. 



1.1 



+0 



m c)-0. 



2 (f \+° 
2\JS)-0 

+0 

-0 



Kfs) 



)±ao(™o) 



-0.2 



(5f 3 )±o: ^o) 



B r 



K+ 



a 



1.6 



K+f 



+o. 

o.9i°; 

1.81°; 

1+0. 
-0. 



g(m c ): 



(m c )18 



2(/s)-0 

i(/s)i° 

2/? ^,+0 
2USJ-0 



+0.6/ K \+0A ( nS ^1+0.0/ w f q\ 
7l%,2]-0.4^1,3)-0. 2 l TO o)(,Io j 



n q+o- 5 ^ "t+o.o/r ^+ - 4 ^^ v 

U - j -0.2( m cJ_o.l(,/sJ_o. 2 ( a l,2)- 



3/ K x+0. 
5l a l,2j-0. 
6/ K \+0. 
6l a l,2j-0. 

0. 
0. 



3 )l§:!h)ffl) 

3 )±o.o(^o)(f„ q ) 
)taS(mo)(fg) 



For i? c — > ao(vr,p) decays, their BRs can be read from the Tables HI and IVl (in unit of 10 6 ), 

Br(B c -> a+7T°) = 6.5±|;| , Br(B c -> agO = 3.5±?;| , (68) 
5r( J B c -> a+p°) = 12.71^ , i?r( J B c -> a°p+) = 10.6±|;| , (69) 

where the various errors as specified have been added in quadrature. One could find the rather 
different decay patterns from these theoretical predictions, i.e., Eqs. ([681 l69|) that Br(B c — > 



ir°) > Br(B c — > a%ir + ) while Br(B c ->■ ag p°) ~ Br(B c — > a[jp + ) within the theoretical 
uncertainties. Because f p (fj) ~ 1.6(1.3) x f^, it is evident that Br(B c — >■ aop) > Br(B c — >■ aovr). 
Based on these pQCD predictions of BRs for i? c — >• ao(vr,p) decays, which are within the reach 
of LHCb experiments [38], it is expected that if the observation or the experimental upper limit 
on the decay modes B c — > aoir(aop) are much smaller than the expectation, this might rule out 
the qq structure for the a,Q. 

On the other hand, the isovector scalar meson ao(1450) has been confirmed to be a conven- 



tional qq meson in lattice calculations 39h43|] recently. Hence, the calculations for the ao(1450) 
channels should be more trustworthy. Our results shown in Tables IIIII and IVIII indicate that 
B c — > ao(1450)7r and B c — > ao(1450)/9 have large branching ratios, of order (5 ~ 20) x 10~ 6 and 
(15 ~ 47) x 10 -6 , respectively. A measurement of them at the predicted level will reinforce the 
qq nature for the do (1450). 

For those B c — > a(P,V) decay modes with oo(1450) as one of the final states, the pQCD 
predictions in Tables iLLTl and IViTl show that for the AS = processes B c — > ao(1450)(7r, 7/ ' , p, u>) 
the BRs in SI are much larger than that in S2, however, for the AS = 1 processes B c — > 
ao(1450)iirW, the BRs in SI are much smaller than that in S2, which will be confronted with 
the ongoing and forthcoming related experiments. It is hoped that the precision measurements 
could help us to determine which scenario is favored by the experiments, then the inner quark 
structure definitely. 

For B c — > a(rj,rj') decays, the numerical results grouped in the Tables HI and Hill indicate the 
small differences between B c — > ar\ and B c — > arj' modes, which is mainly because the relevant 
final state mesons, rf\ contain the same component uu + dd, just with the different coefficients, 
i.e., coscf) and sin^. This pattern is very similar to that of B c — > prj^ decays [io| . 
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For the AS" = 1 B c -> aK^ processes, the pQCD predicted BRs are in the order of 10 , 
which is below the reach of the LHCb experiments (~ 10~ 6 ). From these numerical results as 
displayed in Tables HH HVj EH and ETUI one can find that Br(B c -»• a+K^°) ^ 2 x Br(B c -»• 
a °7^(*)+) although for a meson the vector decay constant f a o = 0, which exhibits clearly that 
the contribution is dominated by the odd Gegenbauer moments in the leading twist distribution 
amplitude of the scalar a meson. This pattern is well consistent with that stressed by the authors 



inRef. 10] 



TABLE III: Same as Table U but for the AS = processes of charmless hadronic B c — > 
(ao(1450),^(1430),/o(1370),/ (1500))(7r,if,?7,7j') decays in SI and S2, respectively. 



AS = 




Decay modes 


BRs (1CT 6 ) 


B c 


-> a (1450)+7r° 


21.0±r 7 (mc)l^(/ S )t^(aJ)i^(S 1 s 3 )i^(mo) (SI) 
6.3l*;*(m c )li;*(/ s )tg;l(aJ)t?:|(flf )3 )18;i(mo) (S2) 


B c 


-> a (1450)°7r+ 


11.9^j(m c )^;I(/ s )«i(aJ)i|?(Sf, 3 )^(mo) (SI) 
4.9^ : 6 3(m c )i^(/5)ir 4 K)ti:2(Sf,3)to:4(^o) (S2) 


B c 


a (1450)+r? 


2.7± ;t(m c )i°;f(/ s )lg; 5 (aI)±^(5^)±^(mo) (SI) 
l-0i O o:|(m c )i° : i(/5)l^(aD±gi(^ 3 )±g : O(m o ) (S2) 


B c 


a (1450) +// 


i.ii83Wi8i(/fi)i85(^)i8:Ws)185(Tiio) (si) 
o.6i8J(fn„)i8:?(/fl)i8:i(^)i8:S(B£s)l85(Tiio) (S2) 


B c 


-> Ko(1430)°K+ 


19.2^(m c )±ti(/ s )l^(a^)^(^ 3 )l^(mo) (SI) 
Q-7±°oi(m c )^(/ s )l°;«(aff 2 )l^(Bf )3 )l°;;(rn ) (S2) 


B c 


-> K°K5(1430)+ 


4.3l oj(m c )l^(/ s )^(a^)i^(Bf )3 )l°; 2 o(rno) (SI) 
9.2l^(m c )t|°(/5)l^(af 2 )iIj( J Bf, 3 )in(mo) (S2) 


B c 


-> fo(1370)vr+ 


3.6^; 2 7 (m c )tr8(/5)irlK)t 1 ^(Bf,3)^^(^o)(fo q ,Sl) 

i-otr 3 (™ c )+o: 2 2 (/5)tr4(«2)t^(sf,3)^:i(^o)(fo q ,s2) 


B c 


-> f O (1500)7T+ 


3.7t^(rn c )tr8(/5)tr1(a2)tij(Bf, 3 )^: 2 3(^o)(fo q ,Sl) 

o.9tr 3 K)-o: 2 2(/5)toi(«2)to:6(sf,3)^:o(^o)(^ 



B. B c -> k(P, V) and B c -> if *(1430)(P, V) decays 

In this type of the considered decays, there are 8 B c — > KqK^*\AS = 0) modes and 16 
B c ->■ K^tt,!]^, p,u,4>)(AS = 1) channels. 

In the AS = processes, we have 4 B c n + K { * ] °, k°K^ + channels in SI and 4 B c 
K^(U30) + K { * )0 , K* (U30)°K^ + decays in both SI and S2, respectively. From the pQCD 
predictions for these considered modes as given in the Tables [U IIIIl PVT and IVIIl one can observe 
that all the BRs are in the range of 10 -6 ~ 10~ 5 within the theoretical errors, which could 
be measured by the near future LHCb and Super-B experiments operated at CERN and KEK, 
respectively. 

Here, it is very interesting to note that for B c — > KqK^ channels Br(B c — > Kq K + ) > 
Br(B c -> K^ + K°) and Br(B c -> K* °K* + ) > Br(B c -> K$ + K*°) in SI, respectively, while the 
situation is quite the contrary for B c — > Kq(1430)K^ decays in S2. One can also find that for 
the B c -> K* (U30)°K^ + decays in both scenarios Br(B c -> K* (U30)° K^+) S1 » Br(B c -> 
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TABLE IV: Same as Table U but for the AS = 1 processes of charmless hadronic B c 
(a (1450), #5(1430), /o(1370),/ (1500))(7r, K, 77,77') decays in SI and S2, respectively. 



AS = 1 
Decay modes 



BRs (10- 



K2)±li(Bls)±°oi(rn ) (SI) 
« 2 )l?:9(^3)±8i(mo) (S2) 

« 2 )^: 6 4(5i s 3)to:i(™o) (si) 

:< 2 )±i:o(£f,3)±o.i(^o) (S2) 

' n\+n.7/ t-..<? N^nn/ » /„.n 



B c -> a (1450) + K° 
B c a (1450)°K+ 



2-3t 
+2 
1 

+0 



6.7^ 
1.2 
3.4 




7 

6 
8 
5 
4 

-1.3 
-0.9 



(fs): 



-0.5 

+ 1.6 (/s)+ 2.3 



-1.4V 

-0.3 (fs)~- 
1.8/ J \- 

-o^Ugj- 



-0.8 



0.6 



:^ni(B?,s)tVo(rno) (SI) 

'^)i^(^)18:8(mo) (S2) 

:4)t°oi(Bl 3 )t o (m ) (SI) 



B c ^K5(1430)+77 



B c ->Kj5(1430)+77' 



5.4+? 
+2 

3.3^ 



6.2^ 



5.1 



+0 



+1 



li(fs) 



+1.5 
-1.6 

-1.3(/s)-0.0 
+0.5/ 7\+0.2 
-0.7\JS)-0.3 
+ 1-2/ J \+0.7 
l.U/SJ-0.3 



ft 



+1.2/7 s+0.4 
1. 3 (/Sj_0.8 
,+0-2/ J -,+0.0 

+0.6/7 -,+0.2 
-0.6\JS)-0A 
,+0.1/7 1+0.0 
-o.iwff-l-o.i 



'i 5 3)±o:o(^o) (S2) 

S.1+0.1/_,_„\ /CM-, 



B c KS(1430)°7r+ 
B c K^(1430)+7r 



6.5 
0. 
3.2 
0.4+° 



+3 
2 
+0 
-0 
+1.1 



+0.3/7 \+0.9 
-0.2US7-Q.4 
+0.6/7 \+0.6 
'— 0.5\J 3 ) -0.3 

'-1.7USJ-2.2 
1+0-1 (f g -1+O.8 

,+0.2/7 \+0.8 
'-0.2 US/ 1 -0.4 
1+0-5/7 \+0.4 
'-0.6USJ-0.4 
,+2-0/7 -,+3.5 
'-1.7USJ-2.3 
+0.2/7 -1+0.8 
-O.lUSJ-0.4 



tt-,+0.9 
"2 J-0.8 

tt-,+0.5 
a 2 J -0.4 

TT-,+0.5 

"2 J-0.4 

TT-,+0.3 

a 2>-0.2 



(*f,3 
(*f.3 



,+o.i 
-0.2 
,+0.0 

'-0.1 
,+0.1 
-0.1 
,+0.0 
'-0.1 



(m ) (SI) 
(m ) (S2) 
(m ) (SI) 
(m ) (S2) 



B c f (1370)K^ 



B c ->■ f (1500)R^ 



0.9 
2.7 



+0 

+0 


7.7+1 
0.7+° 
0.9+° 
2.8+ 

7.9±i" 

>+o. 

>-0 



0.. 



K \+0. 
a l,2)-0 

< 2 ) 



+1 



5/ R S \+0. 
3V-°l,3j-0 
2/ R S \+0. 
ll-°l,3j-0 
6/dS -i+O 



0l £i l,3/'-0 
6/0S -,+0 
5l-°l,3j-0 
K -,+0.5/ n s 1+0 

1,2^-0.: 



if \+0 
a l,2j-0 



5/ dS \+U 
3V-°l,3i-0 
1/dS \+0 
ll-Ol^J-O 
6/0S -,+0 

ol-^i^-o 

,K -,+0.5/oS -i+O 
} l, 2^-0. 5\-°l, 3^-0 



< 2 ) 



+1.' 
2.1 



° (m )(f q ,Sl) 
i(mo)(f q ,S2) 
(m )(fa,Sl) 
° (m )(f s ,S2) 
°(™ )(f q ,Sl) 
i(mo)(f q ,S2) 
4 (m )(fg,Sl) 
g(m )(^,S2) 



Ko(1430)°#W+) S2 , while for the B c -> #5(1430)+ K { *> modes, 5r(5 c -> K*{U30)+K { *> °) S1 < 

(*)0 

5r(B c -> K*(1430)+1T ; ) S2 . It should be stressed that once these predicted BRs and the 
relevant relations could be tested by the experiments in the near future, this could provide 
the great opportunities for us to explore the physical properties of the scalars Kq and the 
corresponding annihilation decay mechanism. 

For the AS = 1 channels B c — > Kq(tt,t]^) and B c — > Kq(p,uj,4>), all the theoretical BRs 
in the pQCD approach are in the range of 10~ 8 ~ 10~ 7 within the theoretical errors except for 
Br(B c -> iT*(1430)(p,w))si ~ 1(T 6 though they are CKM suppressed^ = 0.22), which will 
be confronted by the ongoing and forthcoming relevant experimental measurements. Due to the 
contributions from the same component, i.e., uu, and few differences of the decay constants and 
masses between p° and 00, which result in the similar BRs for B c — > K^p° and B c — > KqU) in 
the considered scenarios. Moreover, we find that the simple relations Br(B c — > Kq°(tt + , /?+)) ~ 
2 x Br(B c — > Kq + (tt° , p )) exists in our pQCD perturbative calculations exactly and Br(B c — > 
K *(1430)(vr,p,o;))si > Br(B c -> if *(1430)(7T,p,u;))s 2 . However, Br(B c -> # *(1430)+</>)si < 
Br{B c —7- Kq (143O)+0)s2j whose pattern agrees well with that obtained by Kim, Li, and Wang 
in Ref. Hj|. 

For B c —7- KQ + (rj,rj') decay modes, based on the pQCD numerical results, we have the 
following remarks: In this sector, both of the components 77^ and r] s in 77 and 77' contribute 
to these channels but with different coefficients even opposite sign. For B c — > k + t]( > de- 
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TABLE V: Same as Table U but for the AS — processes of charmless hadronic B c — > 



i) decays in SI. 


AS' = 
Decay modes 


BRs (10- 6 ) 


B c -> a+p° 
B c -> a°p+ 
B c — > a^w x 10 


12-7^(m c )tii(/ s );^(^)t^(Sf, 3 ) 

10.6l1;S(m c )ll;?(/s)li:§(a§)i?:5(Bf,3) 
9.8t^(m c )t?-j(/ s )t?J(a^)t? : 5 7 (Bf, 3 ) 


B c -> 7t°K*+ 

B c -> K*°K+ 


8.8^i(m c )lJi(/ s )^(af 2 *)l?; Q 7 ( J Bf !3 ) 
4.9±g; 8 7(mc)^:!(/ S )±Ll(aS)i^(^ 3 ) 


B c -> /9+cT X 10 
B c -> /5+fo X 10 


1.6±^(m c )t^(/ s )t^(ag)i^(Sf. 3 )(fo q ) 
0.8^(m c )l°;i(/ g )^(^)l^( J Bf, 3 )(f ^) 



cays, the two parts of contributions make a constructive interference to the branching ratio 
of B c —7- K + r], while a destructive interference to that of B c — > K + rf , which eventually results 
in Br{B c — > n +r n) ~5x Br(B c — > k + t]'). This pattern is very like that of B — )■ K*r\ and K*rf 
decay channels |l|, 0. For B c -> #5(1430)+ r/ ) modes, unlike the £> c — )• K + rj^'\ both of them 
are determined mainly by the factorizable contributions of rj s term, which leads to Br(B c — > 
K£(U30) + i]) ~ Br(S c -»• K^(U30) + i] ) within the theoretical errors in both scenarios. Mean- 
while, it is interesting to note that Br(B c -> K^(1430) + r/('))si < Br{B c -> ^(1430) + r/(')) S2 
while Br(B c — > Kq(1430) + t]) > Br(B c — > Kq(1430) + ?/) in both scenarios, where only the cen- 
tral values are quoted for comparison. Because of the small BRs(< 1CP 6 ) for B c — > Kq + t]( ^ 
decays, all the above theoretical pQCD predictions of the BRs and the physical relations are 
expected to be examined in the forthcoming Super-B experiments. 



C. B c ^f(P, V) and B c -> f(P, V) decays 

As mentioned in the above sections, it is well known that the identification of the structure 
of these neutral scalar mesons / and /' is very difficult, which is a longstanding puzzle not 
yet resolved either by experimentalists or by theorists. Although various scenarios on their 
component have been proposed, by considering the feasibility of factorization approach, we here 
assume these considered scalars to be only qq bound states. 

For the considered 16 B c — > (/, f')(P, V) decays, the numerical pQCD predictions have been 
displayed in the Tables Hl lVIIIl For the / and /', the quarkonia component has been proposed, 
which can be seen in Eqs. ([2]) and (J3j> - For the AS = processes B c — > (f,f')(^ + ,p + ), we use 
the pure qq states /q in / and /' to calculate the BRs in the pQCD approach and obtain the 
numerical results, in which one can find that the BRs of B c — > (tt + , p + )(/o(1370), fo(15QQ))(qq) 
are in the order of 10 -6 within the theoretical errors in both scenarios and within the reach 
of the LHCb experiments (38| . while the BRs of B c — > (ir + , p + )(cr, /o)(qq) are highly below the 
experimental reach of LHCb at CERN. Here, we have assumed that a and /o(1370) have the 
similar decay constant and light-cone distribution amplitudes as /o and /o(1500), respectively. 

As mentioned in Sec. [Hi since the experimental constraints indicate that the mixing angle #0 



between a and /o lies in the range of [25°, 40°] or [140°, 165°] [3l|, then the pQCD predictions 
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TABLE VI: Same as Table U but for the AS = 1 processes of charmless hadronic B c —> 
(a ,K,a,f )(p,K*,uj,4>) decays in SI. 



AS= 1 
Decay modes 



BRs (1(T 7 ) 



B c -> a+K*° 
B c -> agK*+ 



3.5 
1.7 



+0.8 
-0.4 
+0.4 
-0.2 



(m c ) 
(m c ) 



i(fs)±° : 

i(f S n. 



5l M l,2 )-0 



■ 2 (Bl 3 ) 



B c — > K + LU 
B c -> 



1.9 
2.9 



e(m c ) 
e(m c ) 



+o.: 
-o.: 
+o.: 
-o.; 



(fs)tt 



0.6 
0.5 

+O.2c„0N|+O.9 
1.0 



2(«2): 



B c 

B r 



4.5; 

2.3 



t(m c )t° . 



5/1 n+0. 
3ff \+0. 



(sf, 3 ) 
2^-0.4(^1,3) 



0.9 



B, 



B, 



K 



K*+f 



1.6 
2.0 
1.5 
1.9 



i(m c )lg 
l(m c )18 



2/J \+0. 
2USJ-0. 
2/? -i+O. 
2\JS)-0. 

l(/sJ_o.l( a l,2) 



8 (ap + " 
o( a i,2y u 



2V u 'l,2/-0 
5,„if* \+0 



2/J \+0. 
2USJ-0. 



\\ a \:. 



.4(5f, 3 )(f q ) 

i(sf,3)(^) 

i(sf, 3 )(fg) 



of the BRs for B c — > vr + (a, /o) decays with mixing patterns can be read, 
Br(B c -> vr+o-) rj 

Br(B c ^n+f ) « 



(1.9 ~ 2.6) x 10~ 7 for 25° < O < 40° 
(1.9 ~ 3.0) x 10~ 7 for 140° < 6 < 165° ' 



(0.3 ~ 0.8) x 10- 7 for 25° < O < 40° 
(0.1 ~ 0.8) x 10~ 7 for 140° < 9 < 165 c 



(70) 
(71) 



where only the central values are quoted, so are the similar cases in the following text unless 
otherwise stated. Likewise, the pQCD predictions of the BRs for B c — > p + (cr, /o) decays are as 
follows, 



p'a 



Br(B c 
Br(B c -> p+f ) 



(0.9 ~ 1.3) x 10~ 7 for 25° < 9 < 40° 
(0.9 ~ 1.5) x 10~ 7 for 140° < O < 165° 

(0.1 ~ 0.3) x 10~ 7 for 25° < 9 < 40° 
(0.05 ~ 0.3) x 10- 7 for 140° < O < 165° 



(72) 
(73) 



According to Ref. [32], /o(1370) and /o(1500) mixing has the following form, 
/ (1370) = 0.78/ 9 + 0.51/ S , / (1500) = -0.54/ 9 + 0.84/ s , 



(74) 



where we neglect the possible small or tiny scalar glueball components in the present pa- 
per and leave them for future study. Then the pQCD predictions of the BRs for B c — > 
7r + (/o(1370), /o(1500)) decays can be read, 



Br (B c -> 7r+/o(1370)) 
Br(B c -> 7T + / (1500)) 



2.2 x 10- 6 (S1) , 

6.0 x 10~ 7 (S2) ; 

1.1 x 10- 6 (S1) , 
2.7 x 10~ 7 (S2) ; 



(75) 
(76) 
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Likewise, the pQCD predictions of the BRs for B c — > p + (/o (1370), /o (1500)) decays are 

Br(B c -> p+/ (1370)) « { 3 " 7 X ^"^J ' (77) 
V F V " I 1.0 x 10- 6 (S2) ; V ; 



1.8 x 10- 6 (S1) 
5.0 x 10- 7 (S2) . 



Br{B c -> p + / (1500)) « { _ ^_ 7 ^ ' (78) 



For the AS = 1 processes I? c —> K ( -*^ + (f, /') decays, the BRs in the pQCD approach based on 
the pure qq state /q or pure ss one /g of the scalars / and /' are given in the Tables ILT1 II VI IVIl and 
IVIIII One can observe straightforwardly from the tables that all the BRs for B c — > K^ + (f, f) 
channels are in the order of 10~ 8 ~ 10~ 7 except for B c -> if*+/ (1500) in SI though which 
is CKM suppressed. But, if the branching ratio of 10 -6 for B c — > i^T* + /o(1500) decay can be 
detected by the experiments, it is doubtless that the scalar meson /o(1500) is dominated by the 
ss component. When we consider the mixing form for the scalars / and /', the CP-averaged 
BRs for B c —7- (/, /') decays within the pQCD approach have been calculated and shown 

in the Eqs. ([79155]) : 

' (2.0 ~ 2.0) x 10" 7 for 25° < Q < 40° 
(0.5 ~ 1.1) x 10~ 7 for 140° < < 165° 

Br(B c - K+ /0 ) « \ (°- 2 ~ °- 5 ) X 10 " fOT 25 ° < °« < 4 °° ; (80) 
' (0.6 ~ 1.4) x 10~ 7 for 140° < 6 < 165° 



Br(B c — > k + o) « { ; % ; _ : L„ - ( 79 ) 



TABLE VII: Same as Table U but for the AS — processes of charmless hadronic B c 
(a (1450), ^(1430), /o(1370),/ (15QO))(/j, if*, decays in SI and S2, respectively. 



AS = 
Decay modes 



BRs (10" 



B c -> a (1450)+ i o° 
B c a (1450)V 
B c -> a (1450)+w 



+ 16.5 
-12.5 
+ 11.9/ 



(m c )^ 4 2 (/ s ) 



;.3V"2^-7.5 1^1,3 



47.0_ 12 . 5 v'«c;_9.4 

is.s^K)!^^)^;^^)^:?^ 

27.4lg3(m e )jJ»(/ s )^) 
15.5l»;t(m c )H;|(/5)iS:gK) 



6.5+f-(m c )+^(/ s )^-(a^)^-(i?f;3) 

i-ilg:;K)lg:^(/ s )lg:g(^)lg:kgf.3) 



) (si) 

3) (S2) 
) (SI) 
) (S2) 
(SI) 
(S2) 



B c -> K (1430)°K* J 



35.7 
5.4 



-*o 



+ 18 
-12 
+5.5/ 
2.8V 



8/™ \+8.3/ I 

gAjs, 



B c ->■ K K5(1430)+ 



5.o±S:| 
1+2.9 



(m c ) 



3 (m c ) 

" f -)-i:o(/s) 

+1.1 / J \+1.0 



s^'raw ( S1 ) 

,jr"\+5.6/ 



+1 - 4 ^"^ + g:I(o^)r 1 :5(flf, 8 ) (S2) 



Cfe)li:?( 



+2.0 
-1.9 



a h2 )-0.t(-^l,3) 



(SI) 
(S2) 



B c f (1370)p+ 
B c -> f (1500)p+ 



6.ii|?( 

l-7±g;|( 

6-i±i;I( 
1 7+0.5 



3J:5(» 



mc)li5(/5)l?:5( 
m c )ig:I(/5)i^(c 



af)l L |(Bf, 3 )(fo q ,Sl) 

^)ti: 3 (sf, 3 )(f q ,s2) 
^)i?; 2 7 (Bf, 3 )(fo q ,Si) 
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Br(B -> K*+a) « J (3 '° ~ 3 " 5) X 10 ~ 7 f ° r 25 ° < ^° < 4 °° 

\ (0.06 ~ 0.8) x 10~ 7 for 140° < < 165° 

Br(B -+ iT+ M ~ J (0,1 ~ °' 6) X 10_7 for 25 ° < 00 < 4 °° 

1 (2.7 ~ 3.4) x 10- 7 for 140° < 6 < 165° 



Br(B c K+/ (1370)) 
£r(5 c -> if+/ (1500)) 

Br{B c iT+/ Q (1370)) 
£?r(£ c -> iT+/ (1500)) 



1.4 x 10~ 7 (S1) , 
1.8 x 10~ 7 (S2) ; 

7.1 x 10~ 7 (S1) , 
1.3 x 10" 7 (S2) ; 

1.8 x 10~ 7 (S1) , 

6.3 x 10~ 8 (S2) ; 

1.4 x 10~ 6 (S1) , 

5.2 x 10~ 7 (S2) . 



Hence, based on the numerical results shown in the Tables HI HH El and IVH and Eqs. 
2} and (|79M82p . it is evident that the theoretical implications on the components of a and 



TABLE VIII: Same as Table U but for the AS = 1 processes of charmless hadronic B c 
(a (1450),^(1430),/o(1370),/ (1500))(p,^*,w» decays in SI and S2, respectively. 



AS=1 
Decay modes 



B c a (1450)+K*° 
B c a (1450)°K*+ 



B c K5(1430)+o; 



B c -> Kg(1430)^ 



B c -> K5(1430)V 
B c -> K5(1430)+p° 



B c f (1370)K 



*+ 



B c -> f (1500)K 



*+ 



BRs (10~ 7 ) 



2 7 +u - 4 ('rr7 } +u - b 

z - ' -0.6v"cJ-0.5 

7.olr 9 7 K)ii: 7 

l-4lg|K)l^ 
3.5± ;°(m c )+ - 9 



0.8 



(/s 
(/s 
(/s 
(fs 



,+0.1 
'-0.3 
, + 1.2 
'-0.8 
,+0.1 
'-0.1 
,+0.6 
'-0.4 



« 2 *: 

«2 



-0.6 



>±S:!(*?.3 
)±S:l(^f,3 
)iS:S(5? 8 



) (si) 

) (S2) 

) (si) 

) (S2) 



7.1 
1.3 
2.8 
4.3 



3.3 
2.7 
1.2 
0.6 
0.4 
0.6 
+1.7 
1.3 



\+1.5 
cj-l.3 
-,+0.3 
0.2 
0.8 
0.5 



(m c ) 
( 



(fs 
(fs 
(fs 
(fs 



,+0.7 
'-0.8 
,+0.1 
'-0.1 
,+0.3 
'-0.6 
,+0.7 
'-0.4 



(aS?)18:i(i??,3) 

(^)to:s(Bf, 3 ) 
(4)i8:S(^s) 
(4)li:S(fif, 3 ) 



(SI) 
(S2) 
(SI) 
(S2) 



+8.£ 
5.E 
+3.1/ 
1.6k 
2+4.5/ 
■ Z -3.0V 
-1.5/ 
-0.8V 



16.5 
2.9 



1.5 



!(mc)ii: 

-c)^j 
™c)±ij 
" l c)-0.3 



Kfs) 

(fs. 
(fs 
(fs 



+0 

2 

+0.2 
-0.0 
+0.5 
-1.0 
+0.1 
-0.0 



l(a P 2 )t\i(Bl 



P " +2 oi(BL) 

(4)±li mS 



(bL) 



> (si) 

(S2) 
(SI) 
(S2) 



0.3 
0.2 
+0.8 
0.5 

+5 
4 

+2.5 
1.3 
0.2 
0.3 

2.5l°; 7 ( 

3.5±? : g( 



1.0^ 
2.4^ 
14.3^ 
3.2j 

i.r 



(m c )^:Ufs) 

( 



0.3 

0.2 

m c)-0A 

S o(m c )ti 

(mc)i^ 
+0.3 
-0.2 

s+0.5 

l c/-o.5 

N+3. 

c 7-3.3 
-,+0.8 
4 c;_0.7 



(m c ) 



, in 



+0.2 
-0.0 

(/s)tS:S 

Ufs)±l 

+0.5 
-0.4 
+0.0 
-0.1 
+0.4 
-0.3 
+2 

(fs)±£t 



(fs) 
(fs) 
(fs) 

\fsT- 



«2)^:i(Bf.3)(f q J si) 
«2)ti:4(5f,3)(fo q ,s2) 

)(fg,S2) 

)(f q ,si) 

)(f q ,S2) 

3 )(f^si) 

)(fo,S2) 



„( a fS ) -3.3^1,2 

(aS) + -ll(Bf,3 

\ a l,2 J-O.SV-C'l.S 
( a l,2 )-1.6(^1.3 
■ 3 « 2 *)H 4 7(5f, 
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in the light scalar nonet can not be provided by the small pQCD predictions on the short- 
distance contributions of B c — > (ir + , K + , p + , K* + )(a, fo) decays. However, once the large BRs 
above 10~ 6 for AS = processes B c — > (ir + , p + )(/o(1370), /o(1500)) in both scenarios and 
AS* = 1 B c — > A'* + /o(1500) decay in scenario 1 could be measured in the ongoing LHCb or 
the forthcoming Sper-B experiments, they may help determine the components, the ratios of 
quarkonia, and the preferred scenario by the experiments for these two considered scalar /o(1370) 
and /o(1500) mesons, respectively. 

Frankly speaking, for many considered pure annihilation B c decays with BRs of or below 10~ 7 , 
it is still hard to observe them even in LHC due to their tiny decay rates. Their observation 
at LHC, however, would mean a large non-perturbative contribution or a signal for exotic new 
physics beyond the SM. It is worth of stressing that the theoretical predictions in the pQCD 
approach still have large theoretical errors induced by the still large uncertainties of many input 
parameters. Any progress in reducing the error of input parameters, such as the Gegenbauer 
moments of the pseudoscalar or vector mesons distribution amplitudes, Bi of the scalar mesons 
distribution amplitudes and the charm quark mass m c , will help us to improve the precision of 
the pQCD predictions. We do not consider the possible long-distance contributions, such as 
the rescattering effects, although they should be present, and they may be large and affect the 
theoretical predictions. It is beyond the scope of this work and expected to be studied in the 
future work. 

V. SUMMARY 

In summary, we studied the two-body charmless hadronic B c — > SP, SV decays by employing 
the pQCD factorization approach based on the factorization theorem. These considered decay 
channels can occur only via the annihilation diagrams and they will provide an important testing 
ground for the magnitude of the annihilation contributions and implications to the mechanism 
of annihilation decays. Based on the assumption of two-quark structure of the light scalars, 
we make the theoretical predictions on the CP-averaged branching ratios of considered B c — > 
SP, SV channels. In turn, we could obtain the implications on the component and physical 
properties of the light scalar mesons through the experimental measurements on these considered 
charmless hadronic B c decays. Furthermore, these decay modes might also reveal the existence 
of exotic new physics scenario or nonperturbative QCD effects. 

The pQCD predictions for CP-averaged branching ratios are displayed in Tables Ul l Villi From 
our numerical evaluations and phenomenological analysis, we found the following results: 

• The pQCD predictions for the branching ratios vary in the range of 10 -5 to 10 -8 . Many 
decays with a decay rate at 10 -6 or larger could be measured at the LHCb experiment. 

• For B c — > SP, SV decays, the branching ratios of AS = processes are basically larger 
than those of AS = 1 ones. Such differences are mainly induced by the CKM factors 
involved: V u d ~ 1 for the former decays while V us ~ 0.22 for the latter ones. 

• Analogous to B — > K*rj^ decays, we find Br(B c — > k + t]) ~ 5 x Br{B c — > K + rj'). This 
difference can be understood by the destructive and constructive interference between the 
r\ q and rj s contribution to the B c — > K + rj' and B c — > K + rj decay, respectively. 

• For B c — > Aq (1430)7/') channels, the branching ratios for these two decays are similar 
to each other in both scenarios, which is mainly because the factorizable contributions of 
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r] s term play the dominant role and expected to be tested by the forthcoming Super-B 
experiments. 

• If do and k are the qq bound states, the pQCD predicted BRs for B c — > ao(vr,/?) and 
B c — > kK^ decays will be in the range of 10~ 6 ~ 10~ 5 , which are within the reach of the 
LHCb experiments and expected to be measured. 

• For the a (1450) and ^(1430) channels, the BRs for B c -)■ a (1450)(vr, p) and B c -)• 
^(1430)^W modes in the pQCD approach are found to be of order (5 ~ 47) x 10 6 and 
(0.7 ~ 36) x 10 -6 , respectively. A measurement of them at the predicted level will favor 
the structure of qq for the ao(1450) and .fTp (1430) and identify which scenario is preferred. 

• Because only tree operators are involved, the CP-violating asymmetries for these consid- 
ered B c decays are absent naturally. 

• The pQCD predictions still have large theoretical uncertainties, induced by the uncertain- 
ties of input parameters. 

• We here calculated the branching ratios of the pure annihilation B c — > SP, SV decays 
by employing the pQCD approach. We do not consider the possible long-distance con- 
tributions, such as the re-scattering effects, although they may be large and affect the 
theoretical predictions. It is beyond the scope of this work. 
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